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ABSTRACT 

We reformulate the adiabatic inflow-outflow (ADIOS) model for radiatively inefficient 
accretion flows, treating the inflow and outflow zones on an equal footing. For purely 
adiabatic flows (i.e., with no radiative losses), we show that the mass flux in each 
zone must satisfy M oc R n with n = 1, in contrast to previous work in which < 
n < 1 is a free parameter but in rough agreement with numerical simulations. We 
also demonstrate that the resulting two-zone ADIOS models are not dynamically self- 
consistent without the introduction of an energy source close in to the central regions 
of the flow; we identify this with the energy liberated by accretion. We explore the 
parameter space of non-radiative models and show that both powerful winds and 
gentle breezes are possible. When small radiative losses (with fixed efficiency) are 
included, any centrally injected energy flux is radiated away and the system reverts 
to a power-law behavior with n < 1, where n falls in a small range determined by the 
fractional level of radiative losses. We also present an ADIOS model for hypercritical 
(super-Eddington) disk accretion, in which the radiative losses are closely related to 
the flow geometry. We suggest that hyperaccretion can lead to either winds or breezes. 

Key words: accretion, accretion discs — black hole physics — hydrodynamics 



1 INTRODUCTION 

Radiatively inefficient accretion onto black holes can occur 
in the limits of both low and high accretion rates. Although 
black holes can hypothetically absorb the gravitational en- 
ergy liberated during accretion (Begelman 1979), in the re- 
alistic case of a rotating accretion flow the accretion tends 
to be strongly inhibited, with much less gas reaching the 
black hole than is supplied at large distances. The reason 
for this is that the outward transport of angular momentum 
also involves the transport of energy, which causes the gas 
far from the black hole to become unbound (Narayan & Yi 
1994, 1995). 

There is ample observational evidence that black holes 
can be "fussy eaters." The Galactic Center black hole ap- 
pears to be accreting at a rate orders of magnitude below the 
Bondi rate (Melia & Falcke 2001, and references therein), as 
are the black holes in some elliptical galaxies with weak X- 
ray sources (Di Matteo et al. 2000; Mushotzky et al. 2000; 
Baganoff et al. 2003). These systems are believed to be ra- 
diatively inefficient because their densities are low enough to 
inhibit cooling. At the other extreme, cooling is suppressed 
in hypercritical (super-Eddington) accretion flows because 
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radiation is trapped by the large optical depth. Hyperaccret- 
ing sources such as the Galactic binary SS433 (King, Taam 
& Begelman 2000; Begelman, King & Pringle 2006), Galactic 
microquasars during outburst (Mirabel & Rodriguez 1999), 
and possibly ultraluminous X-ray sources (King 2008; Glad- 
stone, Roberts & Done 2009) all appear to produce power- 
ful outflows. Numerical simulations of nonradiative accretion 
onto black holes also show this inhibition, either through the 
development of a zonal structure with inflowing gas near the 
rotational equator, sandwiched between outflowing zones at 
high latitudes, or through convective motions internal to the 
flow (e.g., Stone, Pringle & Begelman 1999; Igumenschev 
& Abramowicz 1999, 2000; Stone & Pringle 2001; Hawley 
& Balbus 2002; McKinney & Gammie 2002; Igumenshchev, 
Narayan & Abramowicz 2003; De Villiers, Hawley & Krolik 
2003; Yuan & Bu 2010). 

Blandford & Begelman (1999, 2004; hereafter BB99 and 
BB04) proposed a simple model for radiatively inefficient 
accretion flows, dubbed the "adiabatic inflow-outflow solu- 
tion" (ADIOS). In their steady state, self-similar formula- 
tion, there are inflowing and outflowing zones with equal 
and opposite mass fluxes that vary with radius according 
to M oc R n with ^ n < 1. In addition, there is pre- 
sumably a small constant mass flux that reaches the black 
hole and is accreted — but this need not be modeled in 



© 0000 RAS 



2 Mitchell C. Begelman 



detail. In the one- dimensional, one-zone models developed 
in BB99 (and presented using a different parameterization 
in BB04), the vertically averaged conservation laws and dy- 
namical constraints — applied only to the inflowing zone — 
impose well-defined relationships among the parameters de- 
scribing the model, such as the specific angular momentum 
of the gas, its effective binding energy (characterized by the 
Bernoulli function), and the mass flux index n. Interestingly, 
n is not determined by these relations, but can take on any 
value between and 1. 

In addition to analyzing some of the two-dimensional 
characteristics of the inflow zone (and showing that the ID 
models provide an excellent first approximation), BB04 con- 
sidered the outflow zone in some detail. In particular, they 
showed that outflowing streamlines can be self-consistently 
patched onto the inflow zone at each radius. By assuming 
that the streamlines are adiabatic, laminar, and inviscid, 
they produced a family of self-similar wind models to repre- 
sent the outflow zone. 

While the BB04 models provide a "proof of concept" for 
ADIOS, they do not necessarily provide an accurate depic- 
tion of even the qualitative features of a realistic zonal flow. 
In particular, the assumption of a wind with adiabatic, in- 
viscid streamlines leads to very different kinds of flows in the 
inflow and outflow zones. While the inflow models depict a 
highly turbulent, well-mixed flow, the outflow models con- 
sist of a well-ordered flow pattern, with streamlines at higher 
latitudes reflecting conditions much closer to the black hole 
than streamlines close to top of the inflow zone. For exam- 
ple, since the terminal velocity along each streamline scales 
with the escape speed at its footpoint, the outflow models 
of BB04 exhibit large shears in the vertical direction, with 
the speed approaching "infinity" near the evacuated cone 
that defines the wind's centrifugal barrier (see, e.g., Fig. 6 
of BB04). 

Numerical simulations that show the generic character- 
istics of an ADIOS (e.g., Stone et al. 1999; Hawley & Balbus 
2002) often seem to exhibit much messier outflows than the 
BB04 models. The outflow region appears to be as chaotic as 
the inflow, with the typical speeds at each radius reflecting 
the local Keplerian velocity. Moreover, in adiabatic simula- 
tions one does not see a range of values for n, but rather 
(where this is measurable) a tendency for the effective value 
of n to be close to 1. 

In this paper, we present a new formulation of ADIOS 
in which the inflow and outflow zones are treated on an 
equal footing. Given the success of the ID inflow model as 
demonstrated in BB04, we adopt a one-dimensional model 
for each zone and mainly consider self-similar flow. We show 
that the introduction of a relatively simple two-zone model 
has a dramatic effect on the character of the solutions: in- 
stead of obtaining solutions for any value of n in the range 
^ n < 1, we find that non-radiative solutions must have 
n = 1. Such flows are regulated by a conserved, outward 
flux of energy through the flow — ignored in the BB99 and 
BB04 models — and in fact require the presence of such 
a flux for self-consistency. The source of this flux must be 
the energy liberated by the residual accretion onto the black 
hole. The inflows associated with these solutions can range 
from marginally bound to nearly Keplerian, depending on 
how large the net energy flux is and how it is partitioned be- 
tween the inflow and outflow zones. Outflow solutions come 



in two varieties: powerful winds with speeds of order the 
Keplerian speed and gentle, subsonic "breezes" driven by 
viscous stresses. 

Flows with self-similar radiative losses have mass flux 
indices n < 1 that depend on the rate of energy loss. These 
flows also depend on the existence of a central energy flux, 
but the net flux in this case is gradually eroded by radiative 
losses. Taking the limit as the radiative loss rate goes to 
zero, we find that these solutions also approach n = 1. The 
inflows in these cases are always relatively tightly bound. 

Finally, we present two-zone ADIOS models for hyper- 
accretion, i.e., super-Eddington accretion, in which the ra- 
diative losses are tied closely to the geometrical and dynam- 
ical structure of the flow. The special form of the energy 
loss rate destroys self-similarity on large scales, but at suffi- 
ciently small radii we can model the radiative effects as small 
corrections to nonradiative flows. By extrapolating these ef- 
fects to scales where they become nonlinear, we develop an 
approximate picture for the flow patterns, luminosities and 
accretion rates of hyperaccreting systems. We suggest that 
there may be two modes of hyperaccretion, corresponding 
to the wind and breeze solutions mentioned above. 

The plan of the paper is as follows. In Section 2 we 
review the one-zone ADIOS model of BB99 and BB04, in- 
troduce the equations for the two-zone model, and show 
that two-zone models without a net energy flux are not self- 
consistent. In Section 3, we explore the parameter space for 
perfectly non-radiative two-zone models and in Section 4 
we consider the effects of (self-similar) radiative losses. In 
section 5 we examine the constraints associated with dissi- 
pative heating. We present our model for hyperaccretion in 
section 6. Finally, we discuss our results and summarize our 
conclusions in Section 7. 



2 ONE- AND TWO-ZONE ADIOS MODELS 

2.1 Review of the one-zone ADIOS model 

We first review the "one-zone" ADIOS model proposed by 
BB99 and BB04, referring readers to Section 2.1 of BB04 
for further details. Adopting the conventions that the local 
accretion rate M and outward viscous couple G = MQ are 
positive (where Q is used in BB04 to represent the viscous 
couple per unit mass), we obtain the following equations 
expressing the conservation of mass, angular momentum and 
energy in a steady state: 



d_ 

dr 



(G-ML) = -(l + r/)L 



dM 
dr 



(i) 



±(GSl-MB) = (f)-l)B*£, (2) 

where L is the specific angular momentum, Q = L/R 2 is the 
angular speed, and 

2 

(3) 



2B? R 

is the Bernoulli function when GM = 1 (where we use a 
Roman typeface for the gravitational constant to avoid con- 
fusion with the viscous couple) and H is the specific enthalpy 
(= 7P/ (7 — l)p for an ideal gas with adiabatic index 7). The 
terms on the right-hand side of equations and |2]) repre- 
sent the rates at which angular momentum and energy are 
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lost to the outflow, per unit radius. For power- law solutions 
with M oc R n , L oc R 1/2 , B oc BT 1 , these forms for the 
angular momentum and energy loss are fully general with 
the appropriate choice of the constants 77 and /3. 

In addition to the conservation laws, the inflow is sub- 
ject to a dynamical equilibrium condition based on the ra- 
dial component of the momentum equation. Given the self- 
similar scalings for pressure and density, P oc i^ n_5 / 2 and 
p oc Jl n ~ 3 / 2 (BB99), and assuming that the inflow is very 
subsonic, we obtain 



i? 2 



R 



+ (5/2 - 



0. 



(4) 



Eliminating P/p in favor of B, one finally obtains a three- 
parameter family of ADIOS models. The independent pa- 
rameters may be taken to be Li? _1//2 , BR and GR~ 3 ^ 2 , or 
n, 77 and f3, or any combination of these. 

BB04 show that solutions can be found for any value of 
n in the interval ^ n < 1. In this interval, the constants of 
integration for the conservation equations may be neglected 
and one obtains simple relations among the various quanti- 
ties. In particular, one can show that the Bernoulli function 
satisfies 

B = GL 1 ~ n (5) 
MR 2 l-n/3 

Since all quantities on the right-hand side of eq. ^ are 
positive with the possible exception of 1 —n/3, we must have 
P > 1/n in order for B to be negative and the inflow to be 
bound. But one can also show that f3 — ► 1 as n — ► 1. The 
sign of B therefore appears to be undefined in this limit. 
Indeed, if one integrates the energy conservation equation 
in the limit n — > 1, one finds that the left-hand side must 
integrate to a constant, whereas the right-hand side diverges 
logarithmically if ft 7^ 0. We conclude that the constant of 
integration — which physically represents a conserved radial 
energy flux through the disk — cannot be neglected in this 
limit. 



2.2 Two-Zone Adiabatic Inflow- Outflow Models 

Let us now generalize the one-dimensional ADIOS formula- 
tion to include the outflow as well as the inflow. In BB04, 
the outflow was treated in a fundamentally different way 
from the inflow. By patching an adiabatic, inviscid stream- 
line onto the inflow at every radius, BB04 showed how the 
outflowing gas could carry away the required amounts of 
angular momentum and energy. However, the resulting flow 
patterns exhibit enormous shear in the poloidal (0) direction 
and large angular gradients of temperature and density. This 
is because each streamline reflects the conditions prevailing 
at the radius where it was launched: streamlines launched 
near a given radius will lie close to the inflow-outflow bound- 
ary and have speeds close to the escape speed at that radius, 
whereas streamlines closer to the axis will have come from 
smaller radii and have correspondingly larger speeds and 
higher specific entropies. 

Simulations, however, suggest that this view of indepen- 
dent, inviscid streamlines may not be accurate. They suggest 
that the outflow is as chaotic as the inflow and is much bet- 
ter described as being "well-mixed." This suggests that it 
might be better to treat the mean properties of the outflow 



at each radius, rather than trying to keep the properties of 
the individual streamlines separate. In this spirit, we treat 
the outflow using the same parametrization as the inflow. 
Denoting the inflowing (outflowing) zones by subscripts 1 
(2), and keeping the sign conventions from Section 2.1 (and 
BB99, BB04) we write down a set of coupled conservation 
equations: 



(Gi -M1L1) = -G12, 
{G 2 + M 2 L 2 ) = Gi 2 , 

(Giili - MiB!) = -Q 12 - S u 
(G 2 Q 2 + M 2 B 2 ) = Q 12 -S 2 , 



(6) 

(?) 
(8) 
(9) 



dr 
d_ 
dr 
d_ 
dr 
d_ 
dr 

where G12 and Q\ 2 represent the fluxes of angular momen- 
tum and energy per unit radius, respectively, transferred 
from zone 1 to zone 2, and Si and S 2 represent net ra- 
diative losses from each zone. By mass conservation, the 
mass fluxes in zones 1 and 2 must be related according to 
Mi(R) = M 2 (R) + Mo 5 where Mo is a constant represent- 
ing the flux of gas that reaches the central object. We can 
therefore drop the subscript on M 2 and write Mi = M+Mo. 

The angular momentum and energy fluxes transferred 
between zones cancel when the pairs of conservation equa- 
tions are summed. We can integrate the summed pairs of 
equations, to obtain 



M(Li 



Gi + G 2 + F L - M Li, 



(10) 



M(B 1 -B 2 ) = GiQi+G 2 n 2 -F E + I (Si+£ 2 )dR-MoB u (11) 

where the constants of integration —Fl and Fe represent 
outward fluxes of angular momentum and energy, respec- 
tively, and we assume that the viscous stresses vanish at the 
innermost radius Rq. Since we are interested in radii far out- 
side the inner boundary, where M ^> M , we can neglect the 
terms containing Mo in both equations. Since ML is an in- 
creasing function of R, we can neglect Fl as well. Since the 
viscous stress is assumed to be directed outward, eq. (10) 
then implies that L\ > L 2 , i.e., that the inflow is spinning 
more rapidly than the outflow at each radius. 

In the one-zone ADIOS model, BB99 and BB04 ne- 
glected Fe as well, which is justifiable for self-similar so- 
lutions with n < 1. In the two-zone case, however, this as- 
sumption yields the curious result that B\ > B 2 , i.e., that 
the inflowing gas is less bound than the outflowing gas. Since 
the motivation for considering ADIOS models is the ten- 
dency of accreting gas to become unbound if it has no way 
to give up energy or angular momentum, this result seems 
suspect. Indeed, we shall show below that it is also incom- 
patible with the condition of dynamical equilibrium. 



2.3 Self-consistency of two-zone models 

To assess the self-consistency of the two-zone models, we 
introduce a generalized version of the dynamical equilib- 
rium condition. We need to include inertial forces in the 
outflow zone, where the speed (y 2 ) may be close to sonic. 
We continue to assume that the inflows are driven by vis- 
cous stresses at (presumably) very subsonic speeds. When 
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mass is being transferred from the inflowing zone to the out- 
flowing zone, the introduction of slowly moving gas emulates 
a drag force in the outflow region. For a power-law model 
with M oc R n , the resulting deceleration is given roughly by 



v 2 



dM 



-nv\ 



2ttRZ 2 OR R ' ^ 

where we ignored the speed of the inflowing gas. The radial 
momentum equation for zone 2 then takes the form 



(1/2 - n)v 2 2 + jj| 



R 



+ (5/2 - n)a 2 = 0, 



(13) 



where we have defined an isothermal sound speed a — 
(P/p) 1//2 . The radial momentum equation for zone 1 takes 
the form 

|| - i + (5/2 - n)ai = 0, (14) 

similar except for the neglect of dynamical terms. As in 
BB04, we ignore the angular correction terms associated 
with the finite scale height of each zone, effectively assuming 
that each zone is geometrically thin. This will lead to errors 
of order (H/R) 2 , where H is the scale height. By compar- 
ing the one-dimensional models to two-dimensional (gyren- 
tropic) models, BB04 showed that the ID models provide an 
excellent approximation, at least for the inflow zones. Since 
the outflow zones sandwich the inflow, we are also implicitly 
assuming that the Hi is substantially larger than H± . 

As in the one-zone models, we can eliminate one of the 
variables in each zone in favor of the Bernoulli function. We 
modify the Bernoulli function for zone 2 to include the radial 
kinetic energy: 



7? - Vl + Ll 



R + j 



7 



2 

a 2 , 



(15) 



with a similar expression (absent the dynamical term) for 
zone 1. If we subtract eq. (13) from eq. (14) and eliminate 
a\ — a\ , we obtain 



Bi — B 2 



(I "») -7(|-n) 



R 2 2( 7 -1)( 



+ 



(f -7(f +») 2 
2(7 " 1) (f " n) " 2 ' 



(16) 



Since the coefficient multiplying L\ 



- L 2 is always negative, 

- L2 must have opposite 
signs if dynamical terms are negligible. This contradicts the 



eq. ( 16 ) shows that B\ — B2 and L\ 



conservation laws, equations (10) and (11), which demand 



that these quantities have the same sign in the zero-flux 
(F E 0) limit. 

In order to change the sign of B\ — B2, the term pro- 
portional to v\ must be sufficiently positive and large. A 
necessary condition is 

ii7) 

which restricts n to values between (for 7 = 5/3) and 1/2 
(for 7 = 1). To check whether consistency is possible even 
for such low values of n, we use the momentum equation to 
eliminate v 2 in favor of a\ 



a 2 . We obtain 



B\ — B2 



71 + 



"it- 2 (I 



+ (7-DU-n) ( ^ 



(18) 



Since n < 1/2, the coefficient of L\ — L\ is still negative, 
implying that the other term must be positive. But in the 
coefficient of a\ — a 2 both the numerator and denominator 
are positive, implying that a 2 < a\. This would mean that 
the outflow must be "colder" than the inflow, which seems 
unlikely given that it lies on top of the inflow, is in vertical 
pressure balance with it, and must have a higher specific 
entropy since it consists of adiabatic gas that has already 
passed through the inflow and has suffered dissipative heat- 
ing in the process. 

Thus, it appears that the constraints of dynamical equi- 
librium in a two-zone ADIOS are incompatible with the con- 
straints of energy and angular momentum conservation, if 
the net energy flux passing through the system is negligible. 
This bolsters our original argument, based on the implausi- 
bility of B\ — B2 being positive on energetic grounds. 

How is this paradox resolved? The only reasonable con- 
clusion seems to be that Fe / and that a two-zone ADIOS 
model requires a substantial flux of energy passing through 
the system from small radii. If radiative losses are negligible, 
this flux is conserved. But in the presence of radiative losses, 
this flux must be large enough to both change the sign of 
B\ — B2 at small radii and compensate for the integrated ra- 
diative losses further out. In the limit of negligible radiative 
losses, this implies that MB is independent of radius and n 
is not a free parameter, but must instead be ~ 1. 



3 NONRADIATIVE TWO-ZONE ADIOS 

We first consider the case in which radiative losses can be 
neglected, and therefore set Si = 82 = 0. Compared to the 
one-zone ADIOS model (BB04), the two-zone model would 
appear to be more restrictive since n — 1 and /3 = 0. How- 
ever, two-zone ADIOS models still form a three-parameter 
family because the energy flux Fe is split into a portion 
that travels through the inflow and a portion that travels 
through the outflow. Defining 



G12 = (l + r?)Li 



dM 
dR ' 



(19) 



we can take 77 to be a third parameter representing the "lever 
arm" associated with angular momentum transfer, with a 
permissible range ^ 77 < 1/2. On physical grounds, we 
expect values of 77 > to apply if organized magnetic fields 
mediate the transfer of angular momentum (BB04). 
Denote the dimensionless energy flux by 



RFe 
M 



(20) 



(where we recall that GM = 1) and let e be the fraction of 
the flux that propagates through the outflow. If we define 
the dimensionless variables 
G „ L 



9 



R 3/2 fl: 



RB, 



21) 



MR 1 / 2 ' R 1/2 
with subscripts 1 or 2 as appropriate, then the integrated 
conservation equations ([6|-(|9| become 



9i 



(1 - 277)^1 



92 



(I+77K1 



(22) 
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gih = bi + (1 - e)fi 



-b 2 + syi. 



(23) 



We therefore have four conservation equations relating the 6 
dimensionless variables gi, <?2, ti, £2, bi, and 62- Note that 
M is not an independent variable, since solutions are invari- 
ant under transformations that keep G/M and Fe/M fixed. 
In addition, we have two dynamical equilibrium conditions, 
which we can cast in terms of the Bernoulli function by elim- 
inating the sound speed a and introducing a dimensionless 
radial speed ui = R 1 ^ 2 V2'. 



bi 

&2 



3-7 
3(7-1) 

3-7 
3(7-1) 



3 + 7 ,2 
2(3-7) 1 

3 + 7 e, 

2(3-7) 2 



57-3 2 

6(7-1) 



(24) 
(25) 



We thus have 6 equations for the 7 variables (including 
u 2 ), which formally define a line through the (/i, £, 77) param- 
eter space. However, we can make the assumption that the 
viscous stress, in either zone, is capable of driving only very 
subsonic radial velocities. Therefore, if g 2 is no n- negligible, 
then u 2 is negligible; conversely, a dynamically significant 
outflow speed implies that g 2 ~ 0, for our purposes. We 
can therefore consider two distinct cases, for each of which 
we have 6 variables and 6 conditions. We refer to these two 
cases as "wind" solutions (g 2 = 0) and "breeze" solutions 
(u 2 = 0), respectively. Given a complete set of parameters 
(/x, e, 77), the two-zone ADIOS model should be completely 
determined in each limit. 



3.1 Wind solutions 

As the wind solution is slightly simpler algebraically, we 
start by analyzing it. Setting g 2 = in the dimensionless 
equations, we immediately obtain 



62 = Sjl 



£ 2 = -(l + r 1 )£ 1 



(26) 



Since b 2 > 0, the outflow is always unbound in the wind 
limit. 

The equations for the inflow region are self-contained, 
i.e., they do not depend on quantities from zone 2, implying 
that the same solution applies for both the wind and breeze 



cases. The Bernoulli function in zone 1 is given by eq. (24), 
but an alternate expression can be derived by eliminating 
gi between the energy and angular momentum equations: 



bl = 3 (1 



(l-e)A*. 



(27) 



Comparing the two expressions, we obtain a simple expres- 
sion for £\ in terms of the input parameters: 



3-7 + 3( 7 -l)(l-e)M 
3 + 7 + 2(7-1)(1-2t,)' 



(28) 



Given £f, we can substitute from eq. (26) into eq. (25) and 
solve for u\. 

We are now able to map out the parameter space for 
wind solutions. The three constraints on physical solutions 
are: 

• £1 < 1, to ensure that the rotation is sub-Keplerian and 
the pressure (sound speed) is positive; 

• bi < 0, to ensure that the inflowing gas is bound; and 




Figure 1. Parameter space of wind models in the e — fi plane 
for 7 = 4/3 and 7 = 5/3, with the angular momentum-transfer 
parameter 77 set to zero. Here and elsewhere, \i is the outward 
energy flux normalized to GMM/R and e is the fraction of that 
flux passing through the outflow zone. For each shaded region, 
the left boundary corresponds to an inflow zone with Keplerian 
rotation (zero-pressure), while the right boundary corresponds to 
marginally bound inflow (61 = 0). The lower boundary of each 
region corresponds to the limit of vanishing outflow speed (u 2 = 
0). Below this boundary, the outflow must be driven by viscous 
stress and the system is in the "breeze" limit. 



• u\ > 0, a self-consistency condition on the wind solu- 
tions. 

These constraints are plotted in Figures [I] and [2] 

Figure [I] shows the parameter space available to wind 
solutions in the e — \i plane, for two values of 7. The angular 
momentum-transfer parameter 77 is set to zero, correspond- 
ing to the purely local transfer of angular momentum from 
inflowing to outflowing gas, i.e., there is no extra moment 
arm in the transfer due, for example, to an organized mag- 
netic field (BB04). The permitted regions are bounded on 
the left and right by the conditions £\ — 1 and b\ = 0, re- 
spectively. This ordering is somewhat counterintuitive, be- 
cause e represents the fraction of net energy flux that passes 
through the outflow. For a fixed total energy flux /x, a larger 
value of e means less energy flowing through zone 1, yet 
the trend is for the binding energy of zone 1 to decrease 
with increasing e. Equally counterintuitive, we note that for 
fixed s, an increase of total energy flux \i drives zone 1 from 
unbound flow toward a thin, cold state. The lower limit of 
each permitted region corresponds to zero outflow speed u 2 , 
clearly a constraint on the validity of the wind model. The 
fact that the allowed parameter space for 7 = 4/3 is more 
restictive than for 7 = 5/3 indicates that it requires more 
energy flowing through the system (and a larger fraction of 
it flowing through the outflow zone) to drive an outflow in 
7 = 4/3 gas. This is not surprising, given that 7 = 4/3 is 
the softer equation of state. 
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Figure 2. Parameter space for wind models with 7 = 4/3 and 
different values of the angular momentum-transfer parameter 77 = 
0, 0.3, 0.5. The more efficiently angular momentum is transferred 
from inflowing to outflowing gas, the lower the central energy flux 
that is required to drive a wind. 



Figure |2] shows the effect of varying 77, with 7 = 4/3. 
As the efficiency of angular momentum transfer increases, 
less energy flux (and a smaller fraction flowing through the 
outflow) is required to drive a wind, in the limiting case 
77 = 1/2, winds can be driven without any net energy flux 
through the system. 



3.2 Breeze solutions 

To study the breeze solutions, we express 62 in terms of £2 
by setting u% — in eq. (25 ). Eliminating gi between the an- 
gular momentum and energy equations for zone 2, we obtain 
a second equation for fr 2 , in terms of £2 and £\ . Eliminating 
62 between these equations, we obtain a quadratic equation 
for £ 2 : 



7 7 - 



•\6 



6(7 - 1) 3 

the solution of which is 



(1 + r?)£i£ 2 + e/* - 



- 7 



3(7-1) 



= 0, 



(29) 



+ 



2(7-1) 
7 7 -3 



±Jq + 



77-3 



3-7 

2( 7 -l)(l + r ? )2 V 3(7-1) 



• (30) 



Since £\ is still given by eq. (28), we can solve for all the 



properties of the breeze solutions in terms of the input pa- 
rameters. 

Breeze solutions are subject to four constraints: 

• i\ < 1 as before, to ensure that the rotation of the in- 
flowing gas is sub-Keplerian and the pressure (sound speed) 
is positive; 



A/1.5 




Figure 3. Parameter space of breeze models in the e — fi plane 
for 7 = 4/3 and 7 = 5/3, with the angular momentum-transfer 
parameter 77 set to zero. For each shaded region, the left bound- 
ary corresponds to an inflow zone with Keplerian rotation (zero- 
pressure), while the right boundary corresponds to marginally 
bound inflow (61 =0) — as in the wind solutions. The upper 
boundary of each region corresponds to the limit of vanishing 
discriminant in eq. ( |30| >. Above this boundary, no breeze solution 
exists. 



• bi < 0, also as before, to ensure that the inflowing gas 
is bound; 

• #2 > 0, the assumption that the stress points outward 
in the outflow, also ensuring that £2 < |(1 + rj)£i < 1; and 

• £2 has real solutions according to eq. (30), which serves 



as a self-consistency condition for the breeze solutions. 



From eq. (30), one can easily see that the condition #2 > is 
well-satisfied when the discriminant in the expression for £2 
vanishes. One can also show that this condition continues to 
be satisfied at all values of the parameters for which a breeze 
solution exists. Therefore, the third condition is superfluous 
and the parameter space of viable breeze solution is deter- 
mined by the conditions on £1, bi, and the existence of a 
real solution for £2 . These constraints are plotted in Figures 
[3] and [4] All breeze solutions have unbound outflow, 62 > 0. 

As before, each solution space occupies a strip between 
the Keplerian limit at the left and the limit of marginally 
bound inflow at the right. The upper boundary marks the 
locus on which the discriminant of the quadratic equation 



(29) vanishes. No physical breeze solution exist for parame- 
ters above this boundary. By comparing Fig. [3] with Fig. ^ 
one can also note that the parameter space for breeze solu- 
tions overlaps with the parameter space for winds; thus, in 
principle one could have either a breeze or a wind model for 
the same set of parameters fi, £,77. Figure [4] shows the effect 
of varying 77, with 7 = 4/3. As in the wind case, a higher ef- 
ficiency of angular momentum transfer decreases the energy 
flux required to drive a wind. 
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1.0 



0.5 



7=4/3 



0.0 




Figure 4. Parameter space for breeze models with 7 = 4/3 and 
different values of the angular momentum-transfer parameter 77 = 
0, 0.3, 0.5. The overlapping bands show the same trend as for 
the wind models, i.e., that a higher angular momentum transfer 
efficiency reduces the energy flux needed to drive a breeze. Note 
that the upper-right threshold for the parameter space is very 
insensitive to the value of 77. 



3.3 Combined parameter space 

We are now in a position to combine our results for the 
breeze and wind models, to show the different possible types 
of solution and how they relate to the input parameters. In 
Figures [5] and [6] we show the combined solution space for 
7 = 4/3 and 7 = 5/3, respectively. The first thing to notice 
is that wind and breeze regions overlap substantially. More- 
over, it is possible that both rapidly rotating (plus sign in 
eq. [30]) and slowly rotating breeze solutions may be pos- 
sible, although the slowly rotating solutions are permitted 
only in a very narrow band of parameter space. (This con- 
straint is based on the assumption that the outflow rotates 
in the same sense as the inflow.) 

These results provide a complete picture of the two- 
zone ADIOS models when losses are neglected. We will now 
consider models that include self-similar radiative losses. 



4 TWO-ZONE ADIOS WITH RADIATIVE 
LOSSES 

Radiative losses are included on the right-hand side of equa- 
tions ^ and ([9]) by specifying the form of Si and 82. 
To guarantee self-similar scaling for these losses, we fol- 
low Narayan & Yi (1994) and define a parameter /, with 
< / < 1, such that —fGdQ/dR equals the net heating 
rate per unit radius. Since the total viscous heating rate per 
unit radius is —GdQ/dR, we have Si — — ( 1 — fi)G\ dQi/dR, 
with a similar expression for zone 2. Note that / = 1 corre- 
sponds to the nonradiative case we have been considering, 
while f —> implies that all locally dissipated energy is ra- 
diated away. We also specify Q12 = (1 — /3)Bi(dM/dR), as 
in BB04. 

We define the dimensionless constants b and £ as in 



3.0 



2.5- 



2.0 



/il.5- 



1.0 



0.5 



y=4/3 



vind only 



wind or breeze 



breeze only 



0.0 



0.0 



0.2 



0.4 



0.6 



1.0 



Figure 5. Combined fi — e parameter space of solutions for 7 = 
4/3, with 77 = 0. Both wind and breeze solutions may exist in the 
dark shaded area. Slowly rotating breeze solutions (i.e., with the 
minus sign in eq. [30] ) may exist only in the narrow strip between 
the wind-only area and the "wind or breeze" region; elsewhere 
only rapidly rotating breeze solutions are allowed. 
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Figure 6. Same as Figure [5] but for 7 = 5/3. 

g and m are functions of R. The dimensionless form of the 
cooling rate then has the form S\R — 3(1 — fi)gi£i/2 in zone 
1, with a similar form for zone 2. 

Under these assumptions, the conservation equations 
admit power-law solutions of the form 



{rh,g u g 2 } = {A,Ci,C 2 }R~ 



(31) 



eq. (21 ), but now define G — gR 3 ^ 2 and M = mR such that 



where A, Ci, C2, and < q < 1 are constants. We may 
choose the overall normalization arbitrarily, so we set A = 1. 
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Substituting eq. (31) into the conservation equations, 
we can obtain solutions for arbitrary choices of 77, fi and 
f 2 , but for the purpose of this section we can simplify the 
algebra considerably by considering the case 77 = with 
f\—fi — f. We then have 



Ci 



and 



3-2q 



C 2 = 



Ci 



(32) 



(33) 



Defining y = 3 — 2q, we obtain 



(2/ - 3/)(*i - t 2 )(ti + yii) + y(y - 3)(6 2 - h) = 0, (34) 

which appears to be a quadratic equation for y (or q) but 
is actually cubic because b\ and b 2 depend linearly on q. To 



see this, we modify the momentum equations (13) and (14), 
replacing n by 1 — q. We then have 



bi = 



b 2 = 



(7-1X3 + 2,?) 



1 



(7-1X3 + 2(2) 



[3 - 7 - 2<z(7 - 1) 

3 + 7 -2g(7-l) £ 2 



[3 - 7 - 2<z(7 - 1) 



3 + 7 -2g(7-l) f 2 



+ 



5 7 - 3 - 2g( 7 - I] 
2( 7 -l)(3 + 2g) 



(35) 



u 2 2 . (36) 



Given i\ and q, we can calculate (5 from the equation 



1-/3: 



1 



t\ 3(1 - /) - 2q 
bi 2(3 - 2q) 



(37) 



which derives from eq. 

These formal solutions to the equations suffer from the 
self-consistency problem described in section 2.3 unless the 
system carries a central energy flux Fe > to partially com- 



pensate for radiative losses. From eq. (11), however, we see 



that the presence of the energy flux would spoil the assumed 
power-law dependence unless 



F E 



(Si+£ 2 )dR. 



(38) 



R 



In such a critically cooled solution, the system "forgets" 
about the energy flux injected at the center and instead 
carries an independent effective energy flux 



f 

J R 



(£i + S 2 )dR, 



(39) 



which allows self-consistency to be maintained. We further 
note that if Fe were smaller than the energy integral in 
eq. ( 38 ) , then no solution would be possible — cooling would 



exceed the total energy supply. Presumably, any radiative 
system will adjust so that the accretion rate yields a cen- 
tral energy flux at least large enough to balance cooling. 
Conversely, if Fe exceeded the energy integral, then radia- 
tive losses would have a minor effect and the nonradiative 
(n = 1) behavior would be preserved. 

Below we consider critically cooled models in the breeze 
and wind limits. 



l.o- 



0.8 



f 

l-f 



7 = 5/3 
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0.2- 



0.0 
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1 -/ 



Figure 7. Allowed values of the exponent q (M oc R 1 ~ q ), nor- 
malized to the radiative fraction 1 — /, as a function of 1 — / 
for 7 = 4/3 and 7 = 5/3,. For each case, the lower bound corre- 
sponds to vanishing angular momentum in the outflow, ^2=0, 
while the upper bound corresponds to vanishing viscous stress, 
G2 = 0. Note that q vanishes oc 1 — / in the limit of negligible 
radiative losses, thus recovering the case n = 1 (M oc R). 



4.1 Breeze solutions 



To extract the breeze solutions, we set u 2 — in eq. (36) 
but make no assumption about g 2 . Substituting for b 2 — b\ 
in eq. (34) and canceling a factor of i\ — £ 2 , we obtain the 



cubic equation 
(y-Zf)(li+yl2)+y(y-S) 



3 + 7-2^(7-1) 
2( 7 -l)(3 + 2g) 



(*i+* 2 ) = 0,(40) 



where we recall that y — 3 — 2q. Breeze solution must satisfy 
the constraints: 

• -1/2 < 61 < 0; 

• £ 2 > 0; and 

• g 2 > 0. 

The last condition implies £ 2 < £±(1 — y~ x ). Since y — 3 = 



—2q < 0, the second term on the left-hand side of eq. (40) is 
negative and we must have y — 3f > 0, which translates to 
q < 3(1 — f)/2. If bi < as well, then eq. (37) implies that 



l-/3> 



1 



(41) 



The allowed range of q(f) is determined by the second 
and third constraints. The actual value of £\ (or b\) is irrele- 
vant, since from eq. (40) we see that q is a function of £ 2 /£i, 



not of either angular momentum separately. The second and 
third constraints imply < £ 2 /£i < 1 — y~ x \ we map this 
region in Fig. [7] 

One can derive these constraints analytically in the limit 
of negligible radiative losses, /—>•!: 



3(7-1) _q_ 
3 7 +l 1 - 



27(7- 



/ 23 7 -3 



(42) 
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As the radiative fraction 1 — / goes to zero, q does as well 
and the solution approaches the case n = 1. 



4.2 Wind solutions 

For the wind case, we take cj2 ~ 0, implying C2 = 0. Com- 



bining equations ( 32 ) and ( 33 ) , we obtain 

2li(l-g) 

£2 = — , 

3 — 2q 



(43) 



which generalizes eq. (26) for £2 in the case 77 = 0. Equation 
( 34 ) can then be written in the form 



h-ti = 3( „ 1 :A~^ 



■[(0-l)(l -?)-«], (44) 



2q(3 - 2g) 4 

where the second version of the equation comes from eq. (|9| 
with cj2 — 0. These equations are identical to eq. (37). 
Wind solutions must satisfy the constraints: 

• -1/2 < 61 < 0; and 

• ul > 0. 

The latter condition places an upper limit on q(f) that is 
identical to the upper boundary for the breeze zone in Fig.[7| 
i.e., corresponding to G2 = 0. The allowed zone for wind 
solutions extends down to q = for all /. The parameter 
space for breezes therefore completely overlaps with that for 
winds, in the q — (1 — /) plane. 



5 DISSIPATION 

Our two-zone models of ADIOS flows have been based on the 
conservation laws of angular momentum and energy, supple- 
mented by dynamical equilibrium conditions in both the in- 
flowing and outflowing zones. In addition, these zones must 
obey the laws of thermodynamics, which describe how en- 
ergy is distributed among various forms. Consideration of 
the thermodynamics of dissipative heating can impose ad- 
ditional constraints on the possible solutions. Since we have 
already found that models are completely determined by the 
choice of three input parameters (/x, e, 77) — plus, under cer- 
tain conditions, an assumption about whether the outflow 
is a wind or a breeze with either fast or slow rotation — it 
is clear that any further restrictions may lead to constraints 
on the possible values of these parameters. For example, it 
may not be possible for the constant energy flux through the 
flow, or its partition between inflowing and outflowing zones, 
to take on any value — the allowed values may instead be 
dictated by thermodynamics. 

We are basically concerned with estimating the rate of 
entropy generation within each zone of the flow. Under self- 
similar conditions with n = 1, the pressure scales with ra- 
dius asPoc R~ 3 / 2 and the density scales as p oc RT 1 / 2 im- 
plying that the specific entropy function scales as P/p 1 oc 
^-(3- 7 )/2 The entropy thus 

increases with decreasing ra- 
dius, as expected for inflowing gas that undergoes viscous 
dissipation (but does not radiate) as it drifts inward. How- 
ever, the self-similar nature of our two-zone model implies 
that the same radial entropy scaling must also apply to the 
outflowing gas — in other words, the specific entropy of the 
outflowing gas must decrease as the gas flows outward. This 
entropy decrease in the outflowing zone indicates that the 



entrainment of low-entropy gas from zone 1 must outweigh 
the entropy increase due to local dissipation. 

These arguments serve as a warning of the care re- 
quired in attempting to estimate dissipative effects. Setting 
s = ln(P/p 7 ), we write the vertically integrated entropy 
equation for zone 1 in the approximate form 



ai „vctei 2nRa 
7 — 1 on 7 



ds. 
pv z —dz ■■ 
oz 



-G 



<Kh dQ 
dR dR' 



(45) 



where the integral is over zone 1 and we have glossed over 
such subtleties as the difference between the vertical aver- 
age of the sound speed used here (recall that a 2 — P/p) 
and those used in the dynamical equilibrium conditions. 
A similar equation applies in zone 2, except that the first 
term on the left-hand side has the opposite sign. For now 
we ignore the dQ/dR term, which represents energy lost 
per unit radius due to conduction or radiation. We set 
2nRpv z w dM/dR, and assume that the second term on 
the left-hand side can be approximated by 



dM 



7- 1 dR 



As, 



(46) 



where As is the difference in specific entropy between the 
top and bottom of the zone under consideration. Presumably 
As > in both zones 1 and 2. Writing ds/dR = -(3-j)/2R 
in both zones, redefining a so it is dimensionless (a 2 — > i?a 2 ), 
and using the dimensionless variables defined earlier, we ob- 
tain 



Asi 



As 2 



3(7 - 1) gih 



3(7- 



1) 92^2 3-7 

o 1 ^ 



(47) 



(48) 



Thus, we see that in order to preserve the self-similar scaling, 
there must be a large entropy jump across zone 2, but there 
need not be one across zone 1. For the sake of making a crude 
estimate, we will assume that the entropy jump across zone 
1 is negligible, and that all the entropy change due to mass 
transfer between zones occurs in zone 2. In this limit, we 
have 



giii 



(1 



"7 



2(1-3) 3(7-1)' 



(49) 



where we have used the equations of energy conservation 
and dynamical equilibrium in zone 1 to express the left- 
hand side in terms of £\. The entropy jump across zone 2, 
which must be approximately equal to S2 — s\ if we assume 
that Asi ~ 0, is likely to be even larger than the estimate 



given in eq. (48), which only includes the Rcj)— component 



of the stress. In addition to the possibility of a significant 
zcj) stress due to the vertical shear in angular velocity, there 
may be large Rz stress. This will certainly be a factor for 
the wind solutions, where we have already had to take into 
account the drag on the outflow due to entrainment of gas 
with little radial velocity (cf. eq. [l3]). 



Solving eq. (l49| for £ 2 , we obtain 



2(3-7) 



(50) 



2(3- 7 ) + 3(7-l)(l-2r ? )- 

In the limit 77 — »> 1/2, angular momentum loss carries off all 
the accretion energy and we are left with a thin, Keplerian 
disk. For all other cases, viscous heating is important and 
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t\ < 1, as expected. Using eq. (24) to express the inflow 



Bernoulli function in terms of £ 1 , we obtain 



h = 



2(3 -j)r, 



(51) 



2(3- 7 ) + 3(7-l)(l-2i/)- 

For the case 77 = 0, bi = and the disk is marginally bound. 
In this limit, the locus of permitted values of \i and e would 
track the right-hand boundary of the solution space in fig- 
ures [5] and [6] For < 77 < 1/2 the disk is bound and the 
locus lies inside the parameter space. For 77 — > 1/2, the locus 
hugs the left-hand boundary of the solution space and the 
Bernoulli function has the Keplerian value b\ — —1/2. 

It proves interesting to generalize our dissipative model 
to the case where Asi is positive, e.g., due to mixing with 



the hotter gas in the outflow zone. In that case, eq. (49) is 
replaced by the inequality 



2n)t\ 



> 



■ 7 



3(7 



(52) 



2(1-*?) 

Since the left-hand side is a monotonically increasing func- 
tion of £\, this means that additional entropy generation 
leads to a colder, more rapidly rotating inflow, and to a 
more negative binding energy. This surprising result could 
explain why simulations of radiatively inefficient flows tend 
to exhibit inflow zones with clearly negative values of the 
Bernoulli function, rather than values hovering close to zero 
(e.g., Stone et al. 1999). 

The effect of radiative losses is to replace the dimen- 



sionless stress g in equation (47) or (48 ) by fg, according to 
the self-similar model of radiative losses adopted in section 
4. For the inflow zone, we obtain 



2(3 - 7) 



(53) 



2(3 -7) + 3/(7 -l)(l-27,)' 

confirming that increasing the losses (decreasing /) makes 
the inflow colder and the rotation faster. The Bernoulli func- 
tion is 

(3 - 7) [1 - /(I - 2r,)] 



h = 



(54) 



2(3 -7) +3/(7-l)(l -2t,)- 

Thus, as expected, an increase in radiative losses makes the 
inflow more bound. On the fi — e parameter space, this moves 
the possible solutions to the left of the zero-energy curve. If 
11 is fixed, this means that e is smaller, corresponding to 
a smaller fraction of the energy flux passing through the 
outflow. 



6 HYPERCRITICAL ACCRETION 

To apply these ideas in a well-defined context, we consider 
the case of hypercritical accretion, in which matter is sup- 
plied at a rate that far exceeds the Eddington limit. Such 
a flow is radiatively inefficient because of the large optical 
depth; indeed, the radiation becomes effectively "trapped" 
in the flow (advected faster than it can escape) at roughly 
the radius within which the liberated power approaches the 
Eddington limit (Begelman 1979). 

Radiative losses are determined using the vertical (z) 
component of the momentum equation. Since radiation from 
zone 1 has to pass through zone 2, we set £\ — in eq. ^ 
and absorb the radiative losses from zone 1 into the overall 
energy loss term Qi 2 . Net radiative losses are thus given 



by £2. In considering the vertical dynamical equilibrium of 
zone 2 we must take account of the outflow speed, which 
can be appreciable. Modeling the vertical velocity profile by 
V2z — V2Z/R oc z/R~ 3 / 2 , we find that the dimensionless 
vertical acceleration is —u^(z/R)/2. There is also a vertical 
drag force, associated with lifting material out of zone 1 and 
mixing it with the outflowing material in zone 2 (cf. eq. |12| ), 
that gives a deceleration —u^z/R). Inserting these terms 
into the momentum equation, we obtain 



H 2 



■ AnGMc ( u\ 

E2 ~^r- 1 + Ty R 



(55) 



where H 2 is the height of the outflow layer, k is the opacity, 
and we include the GM factor for clarity (elsewhere we have 
taken GM = 1). We have doubled the losses from one layer 
to take account of the fact that two radiating outflow layers 
sandwich the inflow zone. To estimate H2, we integrate the 
vertical momentum equation assuming a uniform density. 
This gives H 2 /R = V2a 2 /(l + ^|/2) 1/2 , where a 2 is the 
dimensionless isothermal sound speed on the equator. We 
therefore have 



• 4ttGMc r- 

£2 — V2a 2 

kR 



1 + 



1/2 



(56) 



We define a dimensionless accretion rate rh by M = 
AtvcRtti/k, and assume that matter is supplied at some large 
radius Rout at a rate M out = 47rcR ou trh ou t/ k. If Rout is taken 
to be the radius at which radiation is marginally trapped 
by the inflow, then ra ou t ~ 0(1). If the inner radius R in <C 
-Rout , then it is possible to develop a series solution in powers 
of a small parameter 5 at radii R <C Rout- Importantly, the 
dimensionless coefficients we have used to characterize the 
flows remain constant to first order in 5 — they vary with 
radius only at 0(S 2 ). In contrast, the mass flux and viscous 
stress develop a logarithmic radial behavior at first order in 
(5, which truncates the ADIOS at Rout- 

To see how this behavior arises, we develop solutions to 
first order for zone 1, anticipating the radial dependence of 
rh (which is governed by the form of £2) by writing 



rria 



6\n 1 



rn a R 



3/2 



gia -gia.8 



(57) 
(58) 



where quantities denoted by a subscript "a" are regarded 
as zeroth order and quantities with a subscript "6" are 
first order. We correspondingly define l\ — £± a + £ib and 
bi = bia + bib. Substituting these forms into the conser- 
vation equations, we recover equations ( [22] ) and (23) as 
zeroth-order relations. To first order in the angular momen- 
tum equation, the terms proportional to \n(R/R ou t) give 
9ia — g'ia — 9i 5 implying that Gi(R) = giR 3 ^ 2 rh(R). The 
remaining first-order terms give the correction to the specific 
angular momentum, 



£ib = 



2, I + 77 n 



(59) 



Identifying gi£i a 
we obtain 



-bia = (1 — s)jji as in the nonradiative case, 



(1 - p)b la = 6(1 - s)fi. (60) 
In the presence of radiative losses, the energy transfer rate 
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Q12OCI — /3isno longer zero, but behaves as a first-order 
quantity. 

To proceed further, we need to use the perturbed ver- 
sion of the radial momentum equation. As before, we have 
p = pa 2 , so that p' / p — a 2 (p' / p + 2a' / a) with a' /a = —1/2R 
to first order. Since m oc pvR 2 with v oc a to first order, the 
density must contain the same first-order corrections as m, 
and we have Rp' / p = —1/2 — S. The equilibrium condition 
is then 



1 + 



(§ + «)«! 



0, 



(61) 



which we solve to first order after expanding i\ and taking 
a i ~ a ia + 2ai a aife. Specializing to 7 = 4/3 in the Bernoulli 
function, we obtain 



a la + 2ai a ciib 



f(i- 



5 1 



1-2?/ 



e 2 



bla+i 



3 7-?? 
8 1 - 2r/ 



(62) 



(63) 



A similar procedure can be carried out for zone 2, 
but here we quote only the important result obtained by 
comp arin g terms in the integrated, summed energy equa- 
tion (n). Noting that G 2 (R) oc R 3/2 rh(R) and that F E = 
47rGMcrh a p/ k = rh a pLE, where Le is the Eddington lumi- 
nosity, we find that 



5 = 



\/2< 



Cl2 



VflaP 



1 + 



1/2 



(64) 



The correction term in eq. (57) is negative, implying that 



the ADIOS is quenched at large radii by radiative losses. In 
particular, if we extrapolate the logarithmic correction out 
to radii at which the correction is of order unity (recognizing 
that the approximation is a very crude one in this limit), we 
find that M has a maximum at a radius which we identify 
as Rout, implying 

6- 1 = l + ln( I ^-) = ^. (65) 
V Kin J rn ou t 

For exponentially large values of R ut/Rin, as we expect in 
hyperaccreting systems such as SS433 (where this ratio may 
be ~ 10 3 ), 6 may fall in the range ~ 0.1 — 0.2, justifying our 
approximation scheme. 

To put this on a more quantitative footing, we consider 
solutions for 7 = 4/3, 77 = (as seems reasonable for winds 
driven primarily by the effects of radiation pressure) . Figure 
[8] shows solution curves for eq. ( 65 ) in the e — fi plane for 
different values of m ut, superimposed on the combined pa- 
rameter space for winds and breezes from Fig. [5] All three 
breeze solutions (solid curves) are self-consistent where they 
cross the lower shaded (breeze) region. Wind solutions with 
^lout ^ 0.75 can also be self-consistent, since these cross 
the two upper shaded regions where winds are permitted. 
However, according to our simple model there are no wind 
solutions for m ou t = 1. 

A striking feature of the wind curves is their sensitivity 
to the value of m ou t- As noted earlier, we generally expect 
^lout ~ 0(1), since this is the condition that radiation trap- 
ping be marginally effective at the outer edge of the ADIOS. 
But it makes a qualitative difference whether the value of 
m ut is 1/2 or 1 — the former case permits either breezes 
or winds, whereas only breezes seem to be compatible with 




n 1.5- 



Figure 8. Hyperaccreting ADIOS models in the e — p plane, su- 
perimposed on the combined parameter space for 7 = 4/3 winds 
and breezes from Fig. [5] Solid curves show breeze solutions for 
(top to bottom) rhout = 0.5, 1, 2, while dashed curves show wind 
solutions for the same values of m ou t- Self-consistent solutions are 
possible only where a wind or breeze curve crosses the correspond- 
ing shaded region. Thus, winds are possible only for rh ou t ^ 0.75, 
whereas breezes are possible for all three values of m ou t- The type 
of solution chosen may depend on the outer boundary conditions 
of the flow. 



the latter. Our results suggest that there may be two dis- 
tinct modes of hyperaccretion — breezes or winds — and 
opens the possibility of a system switching between them. 
We will investigate the implications of this possible bimodal 
behavior in a separate paper. 



7 DISCUSSION AND CONCLUSIONS 

We have presented a new version of the adiabatic inflow- 
outflow (ADIOS) model for radiatively inefficient accretion, 
treating the outflow on a similar footing to the inflow. In- 
stead of matching a laminar, adiabatic wind solution to 
a single-zone model for the inflowing region, we treat the 
outflow as a "well-mixed" region characterized by verti- 
cally integrated fluxes of mass, energy and angular momen- 
tum. This two-zone formulation is simpler than the origi- 
nal ADIOS models presented in BB99 and BB04, and re- 
flects the global conservation laws in a more transparent 
way. More importantly, it avoids the large velocity shears 
and entropy gradients predicted as a function of height in 
the original ADIOS model. If the outflows are highly turbu- 
lent, then we argue that a vertically averaged description of 
the outflow may be more realistic. 

This simple modification of the ADIOS model has fun- 
damental implications for the parameter space of physically 
realizable systems. For perfectly adiabatic (nonradiating) 
systems, it implies that the dependence of the mass flow 
on radius is always described by n — 1, i.e., M oc R n = R, 
in contrast with the original formulation in which n is an un- 
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determined parameter with a value between and 1. More- 
over, a two-zone ADIOS model with a well-mixed outflow 
requires a conserved energy flux to propagate through both 
the inflowing and outflowing zones. Absent this external en- 
ergy supply, there is not enough power to propel the outflow, 
while the inflowing gas is unable to relinquish enough energy 
to become bound (in the sense that its Bernoulli function is 
negative). In the original ADIOS model, where this energy 
flux is set to zero, the inflow is able to remain bound only 
because it interacts with the local streamlines of the out- 
flow, which carry much less energy per unit mass than the 
streamlines originating much closer to the central object. 
But in the two-zone model, large amounts of mass injected 
into the outflow at large R mix with large quantities of en- 
ergy liberated much closer in. 

We identify the conserved energy flux with the energy 
released by the flux of matter (Mo) that actually reaches 
the black hole. We characterize this energy flux by two pa- 
rameters, an overall energy efficiency /x defined by the ratio 
of the energy flux to GMM/R (which is independent of ra- 
dius), and a parameter e, describing the fraction of /x that 
propagates through the outflow. Physical models exist only 
for a fixed range within the interval < e < 1 that narrows 
with increasing ji. If ji could reach values ^> 1, the required 
value of e would approach 1, implying that nearly all the 
energy must propagate through the outflow. We stress that 
\i is not an accretion efficiency in the usual sense, because 
it compares the energy liberated by the accreted flux Mo to 
the binding energy of distant matter that never reaches the 
black hole. If the actual accretion efficiency is e, so that the 
energy flux is Fe — eMoc 2 /(l — e), then the accretion rate 



Mo 



e GMM 

~»^r- 



(66) 



A larger value of ji merely means that equality between 
M oc R and Mo = const, occurs farther from the black hole. 
In this case, there would be a wide radial zone in which Mo 
dominates the mass flow. Such a region would resemble a 
unidirectional accretion flow (similar to the early advection- 
dominated accretion flow [ADAF] models of Narayan & 
Yu [1994, 1995]), rather than an inflow-outflow system, 
and would therefore develop a positive Bernoulli function 
(BB99). For this reason, we expect that two-zone ADIOS 
models are restricted to values of \i <^ 2 or 3. 

Because our model is self-similar, we cannot treat 
boundary conditions in detail. However, we suspect that the 
value of fi (and therefore the accretion rate into the black 
hole) and possibly the value of e can adjust in response to 
the energetic demands of the system. A couple of our results 
suggest this. First, when weak radiative losses are introduced 
into the model in a self-similar way, the mass flux index n 
drops below 1 by an amount that depends on the fractional 
energy loss per decade of radius. The effect of this is that Mo 
determines the energy flux through the flow at small radii, 
but that the flow gradually "forgets" about the inner bound- 
ary with increasing R, i.e., as the centrally injected energy is 
radiated away. The case of hypercritical accretion provides 
a more dramatic example, since solutions with a large ratio 
of outer-to-inner radius must have m a ~ 1 + \n(R ou t / Rin) • 
Apparently, matter supplied to a black hole at a highly 
super-Eddington rate will also be accreted at a highly super- 



Eddington rate, in order to supply enough energy to match 
radiative losses in addition to powering the flow. Neverthe- 
less, both ji and e must fall into fairly narrow ranges in order 
for the flows to be self-consistent. (The parameter /3, which 
describes energy transfer from inflow to outflow, must also 
self-regulate to w 1 for n — 1). 

Despite the highly constrained parameter space, our 
solutions display a striking range of behaviors from gen- 
tle breezes, governed by viscous stress, to transsonic winds. 
Such bimodal behavior extends to hypercritical accretion, 
where the dependence of radiative losses on flow geometry 
and dynamics breaks self-similarity close to the trapping ra- 
dius. We note that Shakura & Sunyaev (1973) foresaw the 
possibility of both powerful winds and gentle outflows in 
their pioneering discussion of hypercritical accretion. 

Both breezes and winds are unbound, in the sense that 
the Bernoulli function is positive, but any significant radia- 
tive losses near the outer radii of the ADIOS could quench 
the weaker outflows. It is possible that breeze solutions rep- 
resent large-scale circulations rather than genuine mass loss, 
in which case the total mass of the system would grow with 
time if matter continues to be introduced from outside. 

The strength of the ADIOS approach is that it relies on 
simple conservation laws and dynamical equilibrium condi- 
tions, which impose clear constraints on the space of pos- 
sible flows. Thus it can provide a useful framework against 
which numerical models or observations can be analyzed. 
Beyond the application of conservation laws and dynami- 
cal relations, the ADIOS idea relies on the assumption that 
some physical mechanism exists to bifurcate the available 
energy and to launch an outflow that is spatially distinct 
from the inflowing gas. In the case of the two-zone models 
presented here, it also assumes that the inflow and outflow 
zones are both well-mixed. These assumptions may not be 
valid, or may apply only in certain types of systems (e.g., de- 
pending on boundary conditions, magnetic field structures 
[see, e.g., Beckwith, Hawley & Krolik 2008, 2009], or the rela- 
tive importance of magnetic and radiative stresses) . But the 
two-zone conservation laws and dynamical relations are very 
general and can be adapted to treat alternative assumptions 
about the flow pattern. For example, one can devise a two- 
fluid model that might capture the physics of a radiatively 
inefficient flow dominated by nonlinear convection, with the 
inflow and outflow zones interpenetrating one another. Such 
a model could represent a "convection-dominated accretion 
flow" (CDAF: Narayan, Igumenshchev & Abramowicz 2000; 
Quataert & Gruzinov 2000). We will present such a model, 
and compare its properties to those of the two-zone ADIOS 
model, in a forthcoming paper. 
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